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ABSTRACT 



We present unequivocal evidence for a huge (~ 80 x 100 x 100 Mpc 3 ) super- 
structure at redshift z — 0.27 in the 3-D distribution of radio galaxies from the 
' TONS08 sample, confirming tentative evidence for such a structure from the 7C red- 

f**«> . shift survey (7CRS). A second, newly discovered super-structure is also less securely 

found at redshift 0.35 (of dimensions w 100 x 100 x 100 Mpc 3 ). We present full ob- 
servational details on the TONS08 sample which was constructed to probe structures 
in the redshift range < z < 0.5 by matching NVSS sources with objects in APM cat- 
alogues to obtain a sample of optically bright (E«R ^ 19.83), radio faint (1.4-GHz 
flux density Si. 4 3 mJy) radio galaxies in the same 25 deg area as part-II of the 
7CRS. Out of the total sample size of 84 radio galaxies, at least 25 are associated 
Oh | with the two ~100 Mpc-scale super-structures. We use quasi-linear structure forma- 

q . tion theory to estimate the number of such structures expected in the TONS08 volume 

5-h ' if the canonical value for radio galaxy bias is assumed. Under this assumption, the 

£jj \ structures represent w 4-5<7 peaks in the primordial density field and their expected 

number is low (~ 10~ 2 — 10~ 4 ). Because the TONS08 survey was designed to follow up 
a previously, if tentatively, identified super-structure in the 7CRS, the probability of 
finding two super-structures in TONS08 is uncertain but, assuming that the tentative 
detection of the z=0.35 super-structure is real, must lie between ~ 10~ 4 and ~ 10 
depending on how representative the TONS08 region proves to be. We show that simi- 
lar structures (with similarly low probabilities) are also found in previous radio galaxy 
redshift surveys but their significance has not been fully appreciated because they have 
been traced by very small numbers of radio galaxies. Fortunately, there are several 
plausible explanations (many of which are testable) for these low probabilities in the 
form of potential mechanisms for boosting the bias on large scales. These include: 
the association of radio galaxies with highly biased rich clusters in super-structures, 
enhanced triggering by group/group mergers, and enhanced triggering and/or redshift 
space distortion in collapsing systems as the growth of super-structures moves into the 
non-linear regime. 

Key words: radio continuum: galaxies - galaxies: active - cosmology: observations - 
cosmology: large-scale structure of the Universe 



1 INTRODUCTION 

Super-structures are the largest known structures in the 
Universe. It seems they have evolved from the rare 
>3g) peaks in the initial density fi eld at recombination 
Hnatvk. Lukash fc Novosvadlvilll995t) . They are identified 
in the local Universe as aggregates of clusters of galax- 
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ies, superclusters or large ( ~ 100 Mpc x 100 Mpc) but 
thin (~ 10 Mpc) sheets of galaxies such as the Great 
Wall (iGeller fc Hudmlll98SD . Previous large-scale surveys 
suggest that the power spectrum of galaxy density has 
considerable strength on scales up to at least 100 Mpc 
JPercival et alJ 1200 ll) . Cosmological biasing predicts that 
large-scale dark-matter fluctuations help push small-scale 
fluctuations ove r the density contrast required for gravita- 
tional collapse jKaiserlll983) . meaning that aggregates of 
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rich clusters are expected to trace super-structures in the 
same way that aggregates of galaxies trace rich clusters. 
Study of these structures is difficult due to the huge ar- 
eas of sky needed to survey them. Radio galaxies allow us to 
cover larger area s efficiently because the y are biased trac- 
ers of the mass ^Peacock fc Doddslll994h . The bias factor 
quantifies how well different types of galaxies trace the un- 
derlying mass. Relative bias factors are also used to deter- 
mine how well one popu lation of galaxies traces another. 
IPeacock fc Doddsl £l994) find a relative bias factor of 1.5 
for low redshift radio galaxies with respect to that of op- 
tical galaxies (compared to 3.5 for c lusters). Recent re sults 
from joint 2dFGRS and CMB data llLahav et al.ll2002fr and 
the a nalysis of the bi-spectrum of the 2dFGRS llVerde et alJ 
2002) show that the bias parameter for optically selected 
galaxies is b w 1 with respect to the underlying dark mat- 
ter, which is consistent with no biasing on scales of tens of 
Mpc. There have been suggestions of luminosity- dependent 
biasing in which intrinsically brighter galaxies are more bi- 
ased. This would not be surprising as larger, brighter galax- 
ies may preferentially form in high density regions therefore 
being "born" biased. Radio galaxies are almost exclusively 
found to be associated with massive elliptical galaxies, which 
themselves tend to reside in clusters. It is therefore not sur- 
prising that radio galaxies are biased tracers of mass in the 
local Universe. 

It was only with the advent of radio surveys such a s the 
NRAO VLA Sky Survey (NVSS; ICondon et al"1ll99Sl) and 
the Faint Imag es of the Radio Sky at Twenty -centimeters 
survey (FIRST: llBecker. White fc Helfandlll995l) which now 
cover most of the sky down to sensitivities of a few mjy 
at 1.4 GHz, that we can compile the large area samples 
of faint radio galaxies that we need to study large-scale 
structure. Radio surveys also have the advantage of suffer- 
ing no dust obscuration, unlike their optical counterparts. 
They therefore provide very clean source selection and probe 
to large redshifts ( for NVSS, the median redshift is ~ 1; 
ICondon et aUl 998). Signatures of clustering and large-scale 
structure are now well esta blished in the project ed distribu- 
tion of these sources (e.g. IfBlake fc Wa ll (2002) and refer- 
ences therein) . 

New galaxy redshift sur veys such as the 2d F Galaxy 
Redshift Survey (2dFGRS) (IColless et al.l 1200 ll) and the 
Sloan Digital Sky Survey (SPSS) dStoughton et all l2002h 
which exploit multi-fibre spectrographs to obtain thousands 
of redshifts can be used to match radio positions with galax- 
ies and hence obtain the largest redshift samples of ra- 
dio galaxies with redshifts yet. This allows studies in 3-D 
which have b egun with early subsamples of s urvey data. 
ISadler et all J2002I) and iMagliocchetti et alJ ^2002) have 
matched the 2dFGRS with NVSS and FIRST positions re- 
spectively to measure the Radio Luminosity Function (RLF) 
at low redshift. These studies reveal two distinct populations 
of radio galaxies at faint (5*1.4 ~ 1 mjy) levels: classical ra- 
dio galaxies and quasars powered by Active Galactic Nuclei 
(AGN) and star-forming galaxies in which the radio emis- 
sion arises through synchrotron emission from relativistic 
electrons accelerated by supernovae. An increasing fraction 
of the radio population is identified with starbursts at de- 
creasing flux densities. Due to their lower radio luminosities, 
these objects are typically at lower redshift than the radio 
AGN. Because of its flux-density limit, almost all NVSS ra- 



dio galaxies at z > 0.1 are associated with AGN rather than 
starbursts. 

This paper is concerned with a new study of the 3-D dis- 
tribution of radio galaxies at NVSS flux density levels. It will 
be structured in the following way. In Sec. 2 we describe the 
selection techniques used in the Texas-Oxford NVSS Struc- 
ture (TONS) redshift survey of radio galaxies: this study is 
based on the results of this survey in the TONS08 region 
chosen to overlap with part of the 7CRS in which tenta- 
tive evidence for a su per-structure has already been found 
jRawlings et all2003D . We present the observational data on 
TONS08 in Sec. 3. Sec. 4 describes analysis of the resulting 
redshift distribution including discussion on the significance 
of redshift peaks, methods of assigning an overdensity to 
the distribution and methods for calculating how rare these 
overdensities must be given an assumed bias. We also quan- 
tify the geometry of the super-structure using a principal 
components analysis (PC A). Sec. 5 is a discussion and the 
main points of the paper are summarised in the conclusions 
in Sec. 6. 

Unless otherwise stated, we assume a spatially flat 
ACDM Universe throughout the paper with the following 
values for the cosmological parameters: Hubble constant: 
H = 70 km s _1 Mpc _1 ; h = Ho/100 = 0.7; matter density 
parameter at z=0: Qm(0) = 0.3; vacuum density parameter 
at z=0: fiA(0) = 0.7; rms density variation in 8 Mpc spheres: 
<7 gh -i = 0.94; shape parameter V = 0.21; and index of the 
power spectrum n=l. 

2 SAMPLE SELECTION 

The TONS08 (Texas-Oxford NVSS Structure 08 h region) 
survey is in one of the areas covered by the 7CRS 
i Willott et ai1l2002fr a nd the TexOx-1000 (TOOT) survey 
i Hill fc Rawlingsll2002l) . '. t covers the region: 08 h 10 m 20 s < 
RA < 08' 1 29 m 20 s and 24°10'00" ^ DEC < 29°30'00" 
(J2000). Unlike the low- frequency selected 7CRS and 
TOOT, the TONS08 survey is selected at 1.4 GHz from the 
NVSS. For objects of typical spectral index, it therefore goes 
to fainter radio flux densities than TOOT (which has a 151 
MHz flux density S151 limit of 100 mjy, corresponding to 
5*1.4 ~20 mjy for radio spectral index 1 a ~0.8). In addition, 
TONS08 has an optical magnitude limit (E ^ R ^19.5) im- 
posed on it. TONS08 was specifically designed to seek fur- 
ther evidence for a giant super-structure found in this region 
of sky by analysis of the redshif t distribution of in the 7C- 
II survey jRawlings et alJl2003l) . The redshift distribution 
in the 7C-II survey revealed a prominent spike at z ~ 0.25 
of 5 radio galaxies from the 7CRS. Two further radio- loud 
objects at this redshift were also found: one was just below 
the 151 MHz flux density (Sisi) limit (7C0811+2838), the 
other was just outside the RA limit (NV0840+2805) . 

By going to fainter radio flux densities (Si. 4 ^3 mjy) 
and optical limits (E £s R «19.5), the TONS08 sample was 
optimised for looking at clustering of objects in a region of 
moderate (0 z ^ 0.5) redshifts. Table. Qgives a summary 
of the various radio galaxy samples that are used or referred 
to in this paper. 

1 a is the spectral index for radio sources where the radio flux 
density Si, oc v~ a , where v is the observing frequency 



Two 100 Mpc-scale structures in the 3-D distribution of radio galaxies and their implications 3 



Survey 


TvDe 


radio flux density limits 


optical ]VIa,g limits 


colour cut 


a,r6a. (sr.) 


ref. 


TONS08 


ROz 


Si. 4 >3 mjy 


E <19.83 


B - R >1.8 


0.00688 


This paper 


T0NS12 


ROz 


Si. 4 >3 mjy 


E a R <19.5 


B - R >1.8 


0.00496 


Brand et al. in prep. 


TOOT16w 


ROz 


Si. 4 >30 mjy 


E ss R <19.5 


B - R >1.8 


0.0151 


Brand et al. in prep. 


TOOT08w 


ROz 


Si. 4 >30 mjy 


E ss R <19.5 


B - R >1.8 


0.02786 


Brand et al. in prep. 


TONS08q 


ROz 


Si. 4 >3 mjy 


E S3 R <19.83 


B - R <1.8 


0.00688 


Brand et al. in prep. 


Lacy 


ROz 


Si. 4 >20 mjy 


17< E S3 R, <20.2 


galaxies from APM 


0.0122 


Lacv (2000) 


TOOT 


Rz 


S151 >100 mjy 


NONE 


NONE 




Hill fc Rawlings (20021 


7CII 


Rz 


Sisi >500 mjy 


NONE 


NONE 




Willott et al. (2002) 


Sadler 


ROz 


Si. 4 >2.8 mjy 


14< Bj <19.4 


AGN 1 


0.099 


Sadler et al. (2002) 


2dF 


Oz 


NONE 


Bj <19.4 


NONE 


0.61 


Collcss et al. (2001) 


NVSS 2 


R 


Si. 4 >2.5 mjy 


NONE 


NONE 


10.36 


Condon et al. (1998) 


FIRST 3 


R 


Si. 4 >2 mjy 


NONE 


NONE 


2.61 


Becker et al. (1995) 


Peacock & Nicholson 


Rz 


Si. 4 >500 mjy 


NONE 4 


NONE 


9.3 


Peacock & Nicholson ("19911 



Table 1. Table summarising the various samples used in and/or referred to in this paper. R denotes radio galaxy surveys, O denotes 
surveys with an optical magnitude cut and z deno tes spectroscopic r edshift surveys. Notes: 1. AGN were identified using Principal 
Component Analysis (PCA) applied to the spectra iFolkes et al.|[l99Sf) . 2. NVSS is 50 per cent complete at 2.5 mjy and 9 9 per cent 
complete at 3.5 mjy. 3. FIRST is 95 per cent complete at the given limiting flux densities. 4. iPeacock fe Nicholson! Il99 J) impose a 
redshift limit of 2=0.1. 



The selection criteria for all TONS surveys were based 
on cross-matching positions of objects in radio and optical 
surveys. An initial selection was made of NVSS targets in the 
chosen area of sky with 1.4 GHz flux densities Sia greater 
than 3 mjy (just above t he completeness limit of Sia ~ 2.5 
mJv; ICondon et aljli"99Sl) . Although the FIRST survey has 
better spatial resolution than NVSS, we used NVSS cata- 
logue positions to do our matching. This is because FIRST 
is more likely to resolve out nearby, low surface brightness 
objects hence causing possible incompletenesses within the 
survey. A total of 1148 sources were selected in our survey 
area. The NVSS positions were matched with APM positions 
of objects with i5-band magnitudes less than 19.5. The APM 
survey is a digitised sky survey of Palomar blue (O) and 
red (E) sky survey plates measured by the Automatic Plate 
Meas uring (APM) machine in Cambridge dMcMahon et all 
2002). These magnitudes can be related to the more widely 
used Johnson magnitudes using: 

E = R (1) 
O = B - 0.12(0 - E) (2) 

jMcMahon et al.ll2002l) . 

To ensure a complete, flux-density limited survey we 
adopted the following procedure. We selected any objects 
with APM and NVSS positions with an offset of 20 arc- 
sec from each other. We then plotted radio contours from the 
FIRST survey over optical POSS-II images and overplotted 
the positions of the APM and NVSS objects, identifying real 
identifications by eye. Although time consuming, it is impor- 
tant to do this as both the radio and the optical positions 
are subject to various complex uncertainties. Fig. shows 
some examples of objects where the offset between the APM 
and NVSS positions was greater than 10 arcsec. If we had 
selected radio galaxies by employing a simple cutoff offset 
between radio and optical positions we would have either 
missed real objects or would have contaminated our sample 
with miss-IDs. 

The NVSS radio survey is thought to have rms posi- 
tional uncertainties of < 1 arcsec for relatively strong point 
sources (Sia > 15 mjy) up to about 7 arcsec for the faintest 



detectable sources JCondon et alJll998ft . TONS08.015 and 
TONS08_647 in Fig. show that in practise it can actu- 
ally be more than this. The synchrotron emission detected 
in radio galaxies comes from the extended lobes associ- 
ated with the jets. Fan aroff & Riley class II (FRII) objects 
jFanaroff fc Rilevlll974l) often have two hot-spots where the 
jets terminate. Due to its lower resolution, the NVSS survey 
tends to find a position close to the centre of these double ra- 
dio galaxies. If the lobes are asymmetric, the position can be 
less accurate and more worryingly, two close sources can be 
blended together and missed from the sample. TONS08_263 
in Fig. shows two radio sources with a very small angular 
separation. These have been mis-classified as the two lobes 
of an extended radio galaxy. It is especially important not 
to miss such close pairs of objects in clustering studies. As 
a further check, we obtained FIRST cutouts over a larger 
area for each NVSS position to check if two FIRST sources 
were blended into one NVSS source. This picked up two 
nearby sources which had already been picked up by eye 
and included in our sample. TONS08_858 in Fig.^shows an 
example of how the APM survey can blend two objects to- 
gether. The survey has a scanning resolution of 1 arcsec. In 
practise, objects too close together may be blended together 
making the positions less accurate. The APM catalogues 
are t hought to be Js99 per cent complete dMcMahon et all 
2002). Although positions are accurate in the APM survey, 
the magnitudes can be out by 0.5 mag. The O and E mag- 
nitudes have been modified using corrections derived from 
comparison of the APM magnitudes to GSC-2 magnitudes 
for each POSS-II plate (R. White priv. com.). As our survey 
extends over two POSS-II plates, we modify the magnitude 
cut from our original target of 19.5 to 19.83 (the new limit 
of the most shallow plate). All following analysis uses these 
corrected values. 

The final selection criteria was that O — E Jsl.8. This 
can be related to the more widely used Johnson filter colour, 
B - R by: 

B-R= 0.88(0 - E) (3) 
jMcMahon et al.ll20o3) . This cuts out all the bluer objects 
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08 h 2 1 m 04° 03* 02 s 01' 00 s 2CT59" 58° 57' 08 h 24 m 29° 28 s 27° 26° 25° 24° 23° 22° 

Right Ascension Riaht Ascension 

Center: R.A. 08 21 020 Dec +26 17 18.6 Center: R.A. 08 24 2543 Dec +24 34 31.8 



Figure 1. Four examples of radio galaxies selected in the TONS08 sample. The galaxies are (top left to bottom right): 
TONS08.015,TONS08.263,TONS08.647 & TONS08.858. The small circle denotes the NVSS position and the large circle denotes the 
APM position. TONS08.015 and TONS08.647 show how the NVSS position can be substantially out if the radio flux density is small. 
The figure of TONS08-263 shows how NVSS has treated what is actually two sources as a single extended source due to its lower resolving 
power. The figure of TONS08-858 shows how the APM position can be wrong when it blends two nearby objects. 



which tend to be stars (which by coincidence lie close to the 
line of sight to the radio source), quasars (which generally 
have much higher redshifts than that under interest) and 
starburst galaxies (which tend to be at low redshifts due to 
their low radio luminosity). 

Fig- HI shows the loci of AGN radio galaxies, starburst 
galaxies and quasars in the optical magnitude versus colour 
plane. The template spectral energy distributions (SEDs) 
were obtained from the Gala xy Isochrone Synthesis Sp ectral 
Evolution Library (GISSEL: liruzual fc Ohardotlll993h . The 
SEDs show the paths of galaxies with different absolute mag- 
nitudes. For the evolved el lipticals, we use a B magnitude 
value of M*=-20.9 for L* Jffill fc Lillvlll99lD . The average 
absolute R magnitude of the TONS08 sample is M«-=-22.89. 
This corresponds to ~ 1.5L* and is shown by the solid line. 

It can be seen that the colour cut will divide the pop- 
ulations fairly cleanly. To check that no AGN type galax- 
ies were missed by the colour cut (either unusually blue 
galaxies or galaxies with incorrect colours on APM plates) 
we also obtained spectra for all objects with O — E $CT8 
(the TONS08q sample; Brand et al. in prep.). Three ob- 
jects with O — E -sCl.8 were moved to the TONS08 sample 
(TONS08_422,TONS08_563 and TONS08_669). The spectra 
of TONS08_422 and TONS08_563 show clear characteristics 
of an old elliptical type population as doe s TONS08_669 
which is in the 7CRS ijWillott et al.1 12002|) . Two objects 
with O - E >1.8 were rejected from the TONS08 sample 
and moved into the TONS08q sample (TONS08_980 and 



TONS08_047). TONS08_980 is identified as a nearby star- 
forming galaxy: it has strong S [II] and [O II] emission lines 
and is at very low redshift. TONS08_047 is in the TOOT08 
survey and is identified as a high redshift quasar. These are 
annotated in Fig. [5] 

The final number of objects in the TONS08 sample is 
84. Information about these sources is shown in Table. |5] 
Four of the five 7C-II objects de fining the 7C-II redshift 
spike discussed bv lRawlings et all i2003l) have also been se- 
lected in the TONS08 sample; 7C0818+29 32, very unusually, 
is too optically faint JWillott et al.i r2002). Of the two addi- 
tional objects, 7C0811+2838 is selected (TONS08J279) and 
NV0840+2805 falls outside the survey limits (but is picked 
up by the wider TONS08w survey). The 7C-II objects which 
match the TONS08 super-structure members are shown in 
Table. 03 The TONS08q sample comprises 40 objects, and 
will be the subject of a forthcoming paper. 



3 OBSERVATIONS 
3.1 Optical spectra 

Optical spectra were obtained during the period October 
2000 - February 2002 on the 2.6m Nordic Optical Telescope 
(NOT) using the Andalucia faint object spectrograph, the 
4.2m William Herschel telescope (WHT) using ISIS, the 
2.7m Smith reflec tor at M cDonald with the Imaging grism 
instrument (IGI) faill et al.ll2002bll . and the Hobby-Eberly 
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Figure 2. The O — E colour versus E Magnitude plot for the 
TONS08 and TONS08q samples. TONS08 objects are squares, 
TONSOSq objects are split into the lower redshift starburst galax- 
ies (stars) and quasars (filled stars). The O — £=1.8 colour cut 
which divides the two samples is shown as well as the E=19.83 
mag cut. Objects denoted by an arrow are such because APM O- 
magnitudcs only go as faint as 21.48 and therefore only have their 
minimum colours shown. Overplotted are template SEDs from 
GISSEL (Bruzual & Chariot 1993). To the right of the colour cut 
are the SEDs of an evolved stellar population. The solid line is for 
a galaxy with the average TONS08 absolute i?(f» E) magnitude 
(-22.89). The dashed and dash-dot-dotted lines are L* and 5L* 
respectively. To the left of the colour cut is an SED of a young 
stellar population with the average TONS08q star forming pop- 
ulation absolute ij(f» E) magnitude (-23.94). To see the redshift 
corresponding to each E(xs R) magnitude, see Fig. HI 



telescope (HE T) using the Ma rcario low resolution spectro- 
graph (LRS) Jffill et alJll998tl . The IGI has recently been 
upgraded. All observations from January 2002 were with 
the new volume phase holographic grism. 

Most objects in the TONS08 survey were identified as 
moderate redshift radio galaxies as expected. Spectra for the 
complete sample can be found in the complete paper at: 

:/ /www-astro. physics. ox. ac.uk/~brand/08_paper.ps.gz 
They are sorted by redshift for clarity. In most cases, red- 
shifts were determined from absorption lines (very few ob- 
jects exhibited emission lines). The estimated redshifts are 
presented in Table. |S] These were checked by performing 
a cross correlation between the spectra and both a rest 
frame composite spectrum obtained from the data (see 
Sec. |3]and Fig.|5J and an evolved population GISSEL model 
.Bruzual fc Charlodll99l) . 

Fig. [jj] shows the redshift obtained by eye against the 
averaged cross-correlated redshift for the TONS08 survey. 
There is extremely good agreement between the two differ- 
ent methods of obtaining redshifts for the vast majority of 
objects. Six objects have strong line emission and although 
lines were taken out for the cross correlation, the redshifts 
obtained by cross-correlation showed poor agreement with 
the redshifts by eye. This is not surprising as these objects 
will be of different spectral type and observations would have 
been cut short once it was clear that a redshift could be ob- 



Figure 3. The by-eye redshift versus both cross-correlated red- 
shifts for the TONS08 sample. The three objects with the most 
unclear redshifts are labelled and denoted as stars, objects with 
strong lines are denoted as plus signs and all other objects are 
denoted as squares. 



tained. In all future analysis, we take the redshift obtained 
by eye for these objects. 

There are only three objects for which we don't have 
an unequivocal redshift estimate. TONS08_252B has a very 
poor quality spectrum and only one of the cross-correlation 
redshift estimations matched that obtained by eye. However, 
as this object is very close on the sky to TONS08_252A, it 
is very probable that the redshifts will be similar; we take 
the redshift estimate to be real. TONS08_149 also has a 
poor quality spectrum. The by-eye redshift estimate does 
not agree with those obtained by cross-correlation. These 
redshift estimates are probably the result of flux calibration 
errors which cause a jump in the spectrum where the red 
and blue arms of the WHT ISIS are combined. TONS08_473 
possibly has a lower redshift. However, the cross-correlation 
estimates agree with the higher estimate so we shall assume 
that this is correct. We note here that none of the three 
objects are in either of the super-structures. Consequently, 
if any of these redshifts are incorrect, it could only boost the 
significance of the super-structures. 

We plot the R magnitudes against (by eye) redshifts for 
the TONS08 sample in Fig. The three objects with the 
most inconclusive redshifts lie in approximately the same 
distribution as all other objects in the survey. Also, the red- 
shifts approximately follow the expected underlying redshift 
distribution. Therefore, we take these redshifts to be correct 
in all following analysis. All following analysis is done on 
the redshifts obtained by eye. The only source with a notice- 
ably faint R magnitude for its redshift, is TONS08_1117 (see 
Fig-BJ- This redshift is fairly secure, with the two methods 
of cross-correlation agreeing with the by eye redshift (Ta- 
ble. |2J . Assuming that the APM R magnitude is correct, 
this may be an unusual source. 

A composite spectrum is shown in Fig. To con- 
struct the composite spectrum, the 50 best TONS08 spec- 
tra were trimmed and shifted to their rest-frame. The 
spectra were then normalised using their integrated flux 
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Figure 4. The R magnitude versus redshift plot for the TONS08 
sample. The three objects with the most unclear redshifts are 
labelled and denoted as stars, objects with strong lines are de- 
noted as plus signs and all other objects are denoted as squares. 
Overplotted a re template SEDs of ev olved stellar populations 
from GISSEL iBruzual fc CharlotlliWal) . The dashed line is for a 
galaxy with the average TONS08 absolute R magnitude (-22.89). 
The solid and dotted lines are L* and 5L* respectively. Note that 
the R magnitudes are only accurate to Ri0.3 mag (see Sec. I4.lt . 
TONS08-1117 is marked as it appears to be unusually faint for 
its redshift. This source is discussed in Sec. 131 
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over the range 4000-4500A and then combined by tak- 
ing their median value. We have c ompared this spec- 
trum to a iBruzual fc Charlotl (|l993f) model stellar pop- 
ulation which formed in an instantaneous burst 2.0 Gyr 
ago. The good agreement shows these objects are con- 
sistent with an old, passively evolving stellar population. 
The prominent features used to identify the redshifts 
iSandage. Sandage fc Kristianl Il97ci) of the radio galaxies 
are over-plotted. 

The majority of the TONS08q sample were low radio 
flux density (5*1.4 ^ 10 mjy), low redshift (z ^0.1) starburst 
galaxies or higher flux density, high redshift quasars (see 
Fig.0. Only one object in TONS08q falls between z=0.11 
and z=0.69. Interestingly this object is at 2=0.283, at the 
first redshift spike in the TONS08 sample. This object is an 
X-ray selected QSO (IWei et al]ll999l) . Although interesting, 
we will not include this in any further analysis here. 



3.2 The 3-D distribution 

A representation of the three-dimensional distribution of the 
TONS08 sample is shown in Fig. 03 Clustering of objects can 
be clearly seen at a redshift of z fa 0.27. There are at least 
13 objects tracing a super-structure with a spatial extent of 
at least 80 x 100 x 100 Mpc 3 (the sense is ARAx ADECx 
Az). The size may be limited by the spatial extent of our sky 
coverage in RA and DEC. There is also a second peak in the 
distribution at z «0.35. This structure is approximately 310 
Mpc away from the lower redshift structure (from centre to 
centre measured in co- moving units) and at least 12 objects 
trace a super-structure with a spatial extent of at least 100 



Figure 6. Representation of the three-dimensional distribution 
of the TONS08 sample between 2=0.2 and z=0. 5 (top). Note the 
effect of the selection function in redshift as plotted in Fig. [7J 
Clear groupings of radio galaxies can be seen at z ^0.27 and 
z rjO.35. Also plotted below are two random samples with a 
selection function from the model redshift distribution (Fig. 171. 
TONS08 looks significantly more clustered than either of the ran- 
domly distributed samples. 

x 100 x 100 Mpc 3 (again limited in RA and DEC by the 
sky coverage). 



4 ANALYSIS 

4.1 A model redshift distribution 

The redshift distribution of the TONS08 sample is shown 
in Fig. |7| Overplotted is a model distribution normalised to 
the same sky area and evaluated for the same optical and 
radio selection criteria as TONS08. 

We obtained the model redshift distribution by inte- 
grating a fitted model distribution function of radio galaxies. 
T his was achieved by usi ng the maximum likelihood method 
of iMarshall et all |W83{) to estimate the best fit parame- 
ters for the bivariate (radio and optical) luminosity function 
(BLF). We required a data set over a sufficiently large area 
to reduce the effects of cosmic variance due to large-scale 
structure. The dat a used to constrain this model were from 
ISa.dleretal. I (12002V This cross matches NVSS sources with 
the first 210 fields observed in the 2dFGRS and covers an 
effective area of 325 deg 2 (13 times larger than TONS08). 
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Figure 5. The de-redshifted composite spectrum for the best 50 spectra in the TONS08 sample. The main characteristic features are 
the 4000 A brea k and the nearby Ca K and Ca H lines. Other dominant features are the G and Mg I absorption lines. Also shown is a 
teruzual fc Charlotll993) model stellar population which formed in an instantaneous burst 2.0 Gyr ago (displaced by -1x10 20 Wm 2 A 
for clarity). Note the lack of strong emission lines in this composite spectrum. 
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Figure 7. The redshift distribution of the TONS08 sample with 
the model redshift distribution overplotted (solid line). The ±1<t 
errors on the model are overplotted (dashed lines). We plot the 
distribution using 28 redshift bins between 2=0 and 2=0.5. 



There are two distinct populations in this sample; AGN (60 
per cent) and star- forming galaxies (40 per cent). We used 
only the AGN as the colour cut in TONS08 cuts out the lat- 
ter population. The BLF was modelled with the following 
function: 



p(L,M b ,z) = A(l+z) q — 



(4) 



where p (L, Mr, z) is the source number per unit co-moving 
volume per unit \og w L per unit M b as a function of radio 
luminosity L, optical absolute bj magnitude Mb and redshift, 



z. A, q, 7 and a were the parameters to be fitted. The best 
fit values for the parameters are as follows: logioA=-28.851, 
logi cr=-0.096, 7=1.710, g=1.876. 

The redshift distribution was obtained by integrating 
the BLF over the selection c riteria using a k-correctio n ob- 
tained from template SEDs jBruzual fc Charlotlll993t) . 

Errors on the distribution were found by performing a 
Monte Carlo simulation. As it is likely that other errors will 
also contribute including incompletenesses in both the radio 
and optical surveys at faint flux densities, we included these 
as follows. The NVSS catalogue is « 90 per cent complete 
at Si. 4=3 mjy dCondon et aljll998l) . There is also a 2 per 
cent systematic surf ace density fluctuation across the sky 
jBlake fc Wallll2002l) . This will be a negligible effect com- 
pared to completeness problems at the magnitude limit of 
the APM (and hence TONS08) survey. Before correcting 
for plate-to-plate variations, the A PM magnitudes have a 
global rms uncertainty of 0.5 mag (iMcMahon et alJ 12002). 
This will be reduced by the plate-to-plate corrections. We 
assume an error of 0.3 mag. This is the difference between 
the magnitude corrections of the two POSS-II plates in the 
TONS08 region. We incorporated this into the error on the 
redshift distribution by determining the R magnitude limit 
as a Gaussian distributed value with mean 19.83 and stan- 
dard deviation 0.3. The Monte Carlo method also allowed 
for errors in the derived model parameters. 

The redshift distribution was then calculated for pa- 
rameter combinations whose maximum likelihood values fall 
within ±lcr. This model is the subject of a forthcoming pa- 
per (Brand et al. in prep.). 

The 2dFGRS has a median redshift of 0.1 jFolkes et alJ 

1999). As a consequence, our redshift distribution will be 
less well constrained at redshifts greater than « 0.3. As a 
check, we obtained a model redshift distribution for the sam- 
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Figure 8. The redshift distribution of the iLac^ l2000l> sample 
with the model redshift distribution overplotted. The ±1<t errors 
on the model are overplotted (dashed lines). We plot the distri- 
bution for 16 redshift bins between z=0 and z=0.8. 

pie obtained bv lLacvl (|200(t ) . Fig. |H1 shows that the model 
provides a good fit to his data, even out to high redshift. It 
also shows two redshift spikes in the data although these do 
not appear to be very significant. We consider this sample 
further in Sec. ECHl 

A comparison of the TONS08 and the model redshift 
distribution in Fig. shows that large-scale structure has 
changed the redshift distribution of our sample significantly 
from that expected for a typical area of sky. The model 
predicts 95±11 radio galaxies in this area of sky. Although 
this is more tha n the 84 observed so urces, the difference 
is not significant. iBenn fc Walll lll995t) have already shown 
that large-scale structure can produce ~ 5 per cent vari- 
ations in the source counts over 5x5 deg 2 regions. The 
signatures of large-scale structure in TONS08 seem clear: 
two prominent spikes at z ~0.27 and z ~0.35 with, perhaps, 
void-like regions at lower redshifts, accounting for the low 
overall number of radio galaxies. 

4.2 The significance of clustering in redshift space 

More than one-half of the total sample lie in the regions 
z=0.24-0.28 and 2=0.33-0.38. This is extremely unlikely to 
arise from Poisson fluctuations associated with sampling the 
expected redshift distribution (Fig. |7J . 

We first explored the significance of clustering with a 
simple method using Poisson statistics in each redshift bin to 
calculate the probability that the number of galaxies could 
be greater than the actual number, given the number pre- 
dicted by the model (if the number of galaxies is greater 
than that predicted by the model) or the probability that 
the number of galaxies could be less than the actual number 
given the number predicted by the model (if the number of 
galaxies is less than that predicted by the model). We did 
this analysis excluding the 4 7C-II objects at redshift z w 
0.27 in order to seek evid ence for red shift spikes independent 
of that seen in the 7CRS jRawlines et all2003Tl . Fig.QUshows 
this probability plotted for different binning intervals. It can 



be seen that different binning intervals make little difference 
to the final result. For 20 bins, the probability in the redshift 
interval 2=0.34 to 2=0.36 is 0.002. Because there are 20 in- 
dependent bins, one expects such a low value in 4 per cent of 
a large set of realisations. In the redshift interval 2=0.26 to 
2=0.28, the probability is 0.003. If this was at a random red- 
shift, we would expect a redshift peak of this size or greater 
in 6 per cent of random realisations. However, because this 
redshift pea k occurs in the same place as independent data 
from 7C-II jRawlings et al.ll2003l) and the TONS08 survey 
was designed to find an independent signal in this particular 
bin, we would expect a peak at this same redshift in only 0.3 
per cent of random realisations. There is also a suggestion 
of a void region between redshifts 0.1 and 0.2. Although in 
void regions the small numbers involved prevent this method 
from giving meaningful statistics, a void region within the 
survey would not be unexpected given the number of voids 
in large surveys such as the 2dFGRS and t he filamentary 
distri bution predicted by N-body simulations djenkins et alJ 
1998). 

We performed this analysis for redshift distributions 
corresponding to the +la error on the model redshift distri- 
bution in Fig. [7] For 20 bins, the probability in the redshift 
interval 2=0.34 to 2=0.36 increases to 0.007. One expects 
such a low value in 14 per cent of a large set of realisations. 
In the redshift interval 2=0.26 to 2=0.28, the probability 
increases to 0.006. We would expect a peak at this same 
redshift in only 0.6 per cent of random realisations (given 
again that we are searching for a signal in a specific bin). We 
also note that our error analysis does not allow fully for in- 
adequacies in the parameterised model. Specifically, the true 
uncertainty in the normalisation at 0.3^ z ^0.5 (i.e. beyond 
the redshift at which the model is directly constrained by the 
2dFGRS data) may well be larger than appears to be the 
case from Fig. [7| For these reasons, we feel that evidence 
for the 2=0.35 super-structure is somewhat less secure the 
super-structure at 2=0.27. 

As a further check, we followed the method of 
ISteidel et alJ <|l998f) to estimate the sig nificance of any num - 
ber counts above Poisson expectations. Isteidel et alJ <ll998l) 
smooth their redshift distribution but we prefer to use our 
model distribution. We again took out the 7CRS radio galax- 
ies to check for independent evidence. We considered each 
N neighbouring galaxies and calculated whether their red- 
shifts were closer together than expected from the Erlangian 
distribution. In order to implement this method, we trans- 
formed all redshifts into a coordinate system in which the 
redshift distribution is flat. We could then search for signifi- 
cant clustering above this flat background. We transformed 
our redshifts into a coordinate system t using 

t(z) = f P(z) dz , (5) 
Jo 

where P(z) is the model redshift distribution normalised so 
that J P(z) dz = 1. 

We then calculated the probability that the galaxies 
would be contained in an interval smaller than the observed 
interval t2 

C=p(t <t2) = I X (N1 + N3-1,N2- 1), (6) 

where N2 is the number of galaxies in interval t2, Nl is 
the number of galaxies with redshifts less than the super- 
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Figure 9. The probability that the number of galaxies in each redshift bin could be greater than the actual number given the number 
predicted by the model (if the number of galaxies is greater than that predicted by the model) or the probability that the number of 
galaxies in each redshift bin could be less than the actual number given the number predicted by the model (if the number of galaxies is 
less than that predicted by the model). The results are plotted for 50 bins (solid line), 40 bins (long dashed line), 30 bins (dashed line) 
and 20 bins (dotted line). The figure on the left shows the results with the four 7CRS galaxies at 2 0.27 omitted from the analysis. The 
figure on the right shows the results calculated for the same data but for the model redshift distribution at the +1<t level from Fig. 171 



structure member with the lowest redshift (in the interval 
il), N3 is the number of galaxies with redshifts greater than 
the super-structure member with the highest red shift (in the 
interv al t3), I x is the incomplete beta function jPress et alJ 
1992) and x = (tl + t2)/(tl +t2+ t3). 

We calculated £ for each set of 7V2 of neighbouring 
galaxies where N2 =2,3,4.... To find the significance of each 
candidate super-structure, we created a set of simulated 
samples which had redshifts drawn from the same model red- 
shift distribution. We then compared the £ of the candidate 
super-structure to that of the simulated data. If the candi- 
date £ was smaller than only 1 in 100 simulated datasets 
then we assigned that super-structure a significance of 99 
per cent. 

We found that the two super-structures are both signif- 
icant at the ^99 per cent level. The intervals that maximise 
the significance of the redshift peaks are 2=0.233-0.285 (27 
radio galaxies) and 2=0.33-0.367 (23 radio galaxies). 

4.3 Calculating the radio galaxy overdensity 

Because TONS08 is only ^80 Mpc wide (at redshift 0.27), 
there is an obvious danger of introducing bias by using the 
redshift intervals which maximise the significance of the 
redshift peaks (which correspond to co-moving depths of 
195 Mpc and 135 Mpc) to estimate the radio galaxy over- 
density. To match the data to theoretical overdensities, we 
chose instead to calculate overdensities in 50 Mpc radius 
spheres centred at the peaks of the redshift spikes. In this 
case we find 13 and 12 radio galaxies in the 2=0.27 and 
the 2=0.35 super-structures respectively. Of the 5 7CRS ra- 
dio galaxies in the larger sample defined by 2=0.233-0.285, 
only two (TONS08_526 and TONS08_669) are in the 2=0.27 
super-structure as defined by a sphere of radius 50 Mpc. We 
note that we have taken the centre of the TONS08 field 
for the centre of the 50 Mpc radius sphere. If we choose to 
maximise the number of radio galaxies in a 50 Mpc radius 



sphere, we could move the angular position of the centre of 
the sphere until we found the maximum number possible. 
In this case we find 15 radio galaxies in both the 2=0.27 
and 2=0.35 super-structures. This would give the super- 
structures higher overdensities. 

Note also that we haven't taken into account redshift 
space distortions. If the super-structure is collapsing, this 
will mean there are actually fewer radio galaxies within a 50 
Mpc radius sphere. The effects of this are discussed in more 
detail in Sec. [5] 

We calculate the radio galaxy overdensity using 

TV 

5 g al = 1 (7) 

JVmodcl 

where N is the number of radio galaxies in the super- 
structure and /V mo dci is the number predicted by the model. 
We determined iV mo dci using a Monte-Carlo method. For 
each realisation, we distributed 84 objects at a random RA 
and DEC within the survey limits and at a random red- 
shift weighted by the model redshift distribution. We then 
counted the number of these objects 50 Mpc away from the 
centre of each super-structure and averaged over the num- 
ber of realisations. We predict there should be an average of 
4.04 and 2.72 radio galaxies in the volume of the 2;=0.27 and 
2=0.35 super-structures respectively. The number predicted 
in the 2=0.27 super-structure volume is more than that of 
the 2=0.35 super-structure even though the model redshift 
distribution is higher at 2=0.35. This is because the 50 Mpc 
sphere takes up less of the survey area at high redshift (see 
Fig. [TBI. 

We assume a Poisson error on the number of radio 
galaxies N found within the 50 Mpc sphere. We also cal- 
culate the probability distribution P(7V mo doi)of the expected 
number of radio galaxies jV mo d e i from the errors on the model 
redshift distribution. We determine the joint probability of 
finding this overdensity OD given the data ( Poisson an d 
model distributions) using Bayes theorem (e.g. ISivieJll996t) : 
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Figure 10. The probability of the radio galaxy overdensity given 
the data (incorporating both Poisson and model errors) for the 
TO NS08 2=0.27 (solid line), the TONS08 2=0.35 ( dashe d line) 
the lLacvl <2000l) 2=0.28 (dot-dashed line) and the lLacvl 12Q 
2=0.35 (dotted line) redshift spikes. 
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P{OD\data) oc P(data\OD) x P(OD). 



(8) 



----- 



P(OD) is simply a prior which we set to 1 if <5 ga i > —1 and 
zero otherwise. This prevents JV from being negative. We 
then marginalise over the the expected number iV mo d c i of 
radio galaxies 



P{data\OD) 



P(N\OD,N model )xP(N model )dN model .(9) 



We find an overdensity of 5 ga i ~ ^.^i=n 02 i 
z=0.27 redshift peak and <5 ga i ~ 3.41±q;4| for the z=0.35 
redshift peak. The errors quoted are 68 per cent confidence 
limits. Fig. HOI shows the probability distributions. 



4-3.1 The Lacy sample 

We also performed the above analysis for the lLacvl (I2000T) 
sample. This sample covers nearly twice the area of TONS08 
in a different region of the sky and has a higher flux-density 
limit (S'i.4=20 mjy) (see Table. Q. We plot the 3-D dis- 
tribution of this sample along with two random samples in 
Fig. 1111 Clustering of radio galaxies can be clearly seen at 
z ~0.3 and z ~0.4, corresponding to the peaks in the redshift 
distribution (Fig.[SJl. The random samples do not appear to 
show the same degree of clustering. 

We find six radio galaxies within a 50 Mpc radius sphere 
at 2=0.28 and four radio galaxies within a 50 Mpc radius 
sphere at 2=0.36. In this case, we move the centre of the 
sphere in the RA direction to find the maximum number of 
radio galaxies within the sphere. This is because the Lacy 
survey is larger in extent in this direction (« 300 Mpc at 
2=0.36) than the TONS08 survey. We are also confident in 
this case that we are tracing the entire super-structures. We 
predict there should be an average of 0.7 and 0.67 radio 
galaxies in the volume of the 2=0.28 and 2=0.36 spheres 
respectively (this is less than obtained for TONS08 due to 
the higher radio flux density limit of the Lacy sample). 
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Figure 11. Representation of the three-dimensional distribu- 
tion of the Lacy sample between 2=0.2 and 2=0.5 (top). Note 
the effect of the selection function in redshift as plotted in Fig.lsl 
Also plotted below are two random samples with a selection func- 
tion from the model redshift distribution (Fig. |gj . Similarly to 
TONS08, this sample looks significantly different to the random 
samples 



We calculate the probability distribution of the overden- 
sity OD given the data for the two redshift spikes in the Lacy 
sample (Fig. llOH . We find an overdensity of <5 ga i ~ 7.57± 1 ;i3 
for the 2=0.28 redshift peak and <5 ga i w 4.97±l\n for the 
z=0.36 redshift peak. The errors quoted are 68 per cent 



confidence limits. We note here (Fig. 1101 that the errors are 
much larger. There is a much higher probability that these 
spikes correspond to lower overdensities than they appear 
which is why, perhaps, not much attention has been drawn 
to them previously. 



4.4 Finding super-structures from radio galaxy 
optical surveys 

Fig. ll2l shows the distribution of apparent magnitudes in the 
TONS08 sample. A "background" level is also plotted from 
a sample selected in the same way over a much (~ 13-times) 
larger area of sky (but normalised down to the same sky 
area). There are no very significant spikes in the magnitude 
distribution. This is because objects at a given redshift have 
a large spread in absolute magnitude (Fig. |1J and therefore 
any structures will become smeared out. This shows that 
redshift surveys are required to pick up significant super- 
structures. 
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Figure 12. Histogram of the apparent magnitudes of the 
TONS08 sample (solid line) against a background level (dotted 
line) obtained from a similar sample in a larger ~ 13-times (but 
then normalised) area of sky. The shaded area shows the 13 ob- 
jects that are members of the 2=0.27 super-structure. 



4.5 Angular clustering 

To test for any angular clustering we applied two- 
dim ensional Ko lmogorov-Smirnov (K-S) tests to our sam- 
ple (iPress et al 11993) . We compared the 2-D distribution of 
radio galaxies in significant super-structures (Sec. 14.21 with 
the distribution of samples randomly distributed over the 
selection volume. To obtain a randomly distributed sample 
for the larger super-structure samples, we distributed the 
sample randomly over the sky area of TONS08. For the sam- 
ples that are defined to be within a 50-Mpc-radius sphere, 
we did this by randomly distributing the sample within the 
volume of the 50 Mpc radius sphere and de-projected this 
back down to a 2-D angular distribution. The Kolmogorov- 
Smirnov D is defined as the maximum value of the absolute 
difference between the two cumulative distributions. This is 
generalised to 2-D by considering the fraction of data points 
in ea ch of four quadrants around a given point dPress et alJ 
Il992l) and taking the maximum difference to define D. 

We compare the real distribution with a random distri- 
bution to find the value of D. To calculate the probability 
that the real distribution is picked from a different distribu- 
tion to that of a random distribution, we also calculate D 
for two random distributions. We do this 10000 times and 
calculate the fraction of random values of D which are larger 
than that of the real D to obtain the probability of the ran- 
dom D being larger than the real D (see Table. If this 
probability is small, the distribution is significantly different 
to a random distribution. 

We find very weak evidence for angular clustering in 
both the super-structures but not within the entire sample. 
The 2-D distributions are shown in Fig. 1131 Evidence for 
angular clustering, at least in the restricted set of super- 
structure members, include an apparent northern edge at 
DEC ~28° in the 2=0.27 super-structure and an apparent 
edge at DEC ~26° in the 2=0.35 super-structure. Edges to 
the distribution can also be perceived in Fig. [(J The symbols 
representing radio galaxies are scaled with redshift to check 



0.27 

0.233-0.285 
0.35 

0.330-0.367 
0.0-0.5 



13 
27 
12 
2.3 

84 



0.38 
0.24 
0.40 
0.31 



0.43 
0.66 
0.35 
0.37 



0.16 0.58 



Table 2. Results of two-dimensional K-S tests for the spatial 
distribution of the significant super-structures (for both the larger 
sample and the sample which includes all radio galaxies within a 
50 Mpc sphere) and the whole TONS08 sample. D is defined as 
the maximum value of the absolute difference between the data 
and a random cumulative distribution. The Probability is that of 
a random-random value of D being greater than the calculated 
value of D. 



there is no systematic velocity shift across the fields. No 
such effect is statistically significant but we note hints that 
the lower redshift objects in both superstructures are more 
tightly clustered than the whole sample. We have found very 
tentative evidence for angular clustering in the samples cor- 
responding to the 50 Mpc spheres. This suggests that both 
these super-structures are real. 



4.6 The geometry of the super-structures 

The spatial extents of our 2=0.27 and 2=0.35 super- 
structures are at least « 80 x 100 x 100 Mpc 3 and 100 
x 100 x 100 Mpc 3 respectively (A RA xA DEC xAz 
in co- moving coordinates). These are limits due to select- 
ing all radio galaxies within 50 Mpc of the centre of each 
super-st ructure. We used Pri ncipal Component Analysis 
(PCA) iMurtaen fc Hecklll987fl to quantify whether the ra- 
dio galaxies trace a spherical or more filamentary distribu- 
tion. PCA seeks the axis from which all the points have the 
smallest variance. The direction of this axis forms the first 
principal component (the first eigenvector). The next prin- 
cipal component is found by minimising the variance in a 
direction orthogonal to the first principal component. This 
process is continued to form three orthogonal axes. We con- 
vert the spatial and redshift coordinates to co-moving dis- 
tances from the centre of the super-structure, making the 
parameters in correct proportion to each other. The spatial 
distributions of radio galaxies in the two super-structures 
and their principal components are shown in Fig. 1141 The 
eigenvalues are the variances of each principal component. 
We can therefore find the preferred direction of the spa- 
tial distribution (the first principal component) and quan- 
tify how much of the variance this accounts for (the eigen- 
values). If the eigenvalues are all fairly equal in size this 
indicates a fairly spherical distribution. Conversely, if one 
eigenvalue is much larger than the others, this indicates a 
filamentary structure in the direction of the corresponding 
principal component. The results in Table. [3] show a slight 
elongation for the 2=0.27 super-structure. The direction of 
the eigenvectors shown in Fig. 1141 shows that this is mainly 
in the 2-direction. 

Although we are performing the analysis on 50 Mpc ra- 
dius spheres, there may be a bias introduced in the preferred 
direction of the PCA analysis. This is because the area of 
the survey corresponds to only about 80 Mpc in the RA 
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Figure 13. The angular distribution of the 2=0.27 (left) and the 2=0.35 (right) super-structure members. The TONS08 survey is 
delimited by the axes surrounding each plot. Note, however, the members of the super-structures (filled circles) are restricted to a 
50-Mpc-radius sphere whose angular projection is overplotted. We also plot the larger number of candidate super-structure members as 
described in Sec. l4.2j Radio galaxies are scaled linearly according to their redshift (lower redshift objects are larger). The redshifts in the 
50 Mpc spheres lie between 2=0.26 and 2=0.281 for the 2=0.27 super-structure (left) and between 2=0.344 and 2=0.366 for the 2=0.35 
super-structure (right). 



direction for the 2=0.27 super-structure. This effect will be 
even more pronounced for the larger samples. To take the 
survey volume into account and calculate the significance of 
any elongation, we performed a Monte-Carlo analysis. 

We calculated 100000 random realisations of the same 
number of objects within a 50 Mpc radius sphere (and within 
the redshift limits for the larger samples) for each super- 
structure. For each principal component, we found the per- 
centage of these realisations whose eigenvalues were as far or 
greater from the mean random eigenvalue as the real eigen- 
value. The results are shown in Table. [3] For the first prin- 
cipal component of the z=0.27 super-structure, only 2.5 per 
cent of the random realisations had eigenvalues as large as, 
or larger than, that measured. The geometry of the survey 
cannot, therefore, be responsible for this high eigenvalue. 

We have not, however, corrected for redshift distortion 
effects. If the super-structure is in-falling, this could make 
the super-structure appear less elongated in the 2-direction. 
This could mean that the real elongation is even more pro- 
nounced. We neglect redshift space distortions here (as they 
will only make the results more extreme) but see Sec. 14.71 
for a more thorough treatment. 

4.7 Implications of the radio galaxy overdensity 

We have established that the 2=0.27 super-structure is sig- 
nificant and the z=0.35 super-structure probably is. We now 
attempt to calculate the probability of the TONS08 survey 
intercepting such structures given the cosmological param- 
eters and different degrees of radio galaxy bias. 

We define the bias factor b as the ratio of the radio 
galaxy overdensity Sp g /p s to the underlying matter over- 
density 8p m /p m 



Redshift 






PCI 


PC2 


PC3 


0.27 


Prop. 




0.637 


0.191 


0.172 




random 


Prop 


0.369 


0.370 


0.260 




MC 




2.5% 


20.7% 


39.7% 


0.233-0.285 


Prop. 




0.669 


0.196 


0.135 




random 


Prop 


0.717 


0.182 


0.101 




MC 




8.6% 


40.5% 


11.1% 


0.35 


Prop. 




0.444 


0.370 


0.186 




random 


Prop 


0.308 


0.345 


0.347 




MC 




80.8% 


57.1% 


20.9% 


0.330-0.367 


Prop. 




0.474 


0.304 


0.222 




random 


Prop 


0.413 


0.386 


0.200 




MC 




19.3% 


41.1% 


21.0% 



8p u 



(10) 



Table 3. The proportion of the total variance (Prop.) ac- 
counted for by the 1st (PCI), 2nd (PC2) and 3rd (PC3) prin- 
cipal components for the 2=0.27 and the 2=0.35 super-structures 
as well as the larger samples. Also shown is the proportion of 
the total variance for random realisations of the sample (ran- 
dom Prop.), and the percentage of random realisations with the 
measured value or further from the mean (MC). The coefficients 
of the principal components are: PC1=[-0.17,0.29,0.94],PC2=[- 
0.97,0.12,-0.21] and PC3=[0.17,0.95,-0.26] for the 2=0.27 super- 
structure and PCl=[0.64,-0.24,-0.73],PC2=[0.55,-0. 53,0.65] and 
PC3=[0.54,0.82,0.20] for the 2=0.35 super-structure 



We assume the local value for the radio galaxy bias fac- 
tor is 6=1.65±0.15. This is determined by multiplying the 
relative bias factor of radio gal axies to optical galaxies 
(1.9/1.3=1.5±0.09 calculated bv iPeacock fc Doddsi ll994T) 
with the bias factor of optical g alaxies relative to the 
dark matter (1.1±0.08 calculated bv lLahav et aljl2002f) and 
adding the errors in quadrature. 
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Figure 14. The spatial distribution of the z=0.27 super-structure (top) and the z=0.35 super-structure (bottom). The spheres repre- 
senting the radio galaxies are scaled logarithmically with their radio flux densities. The solid line shows the first principal component. 
The dashed and dotted lines show the second and third principal components respectively (see Table. [3] for coefficients). The figures on 
the left show the radio galaxies defined as being within 50 Mpc of the centre of the super-structure. The central figures additionally show 
the larger samples (white spheres). The figures on the right show both the larger (white spheres) and smaller samples (black spheres) 
with the principal axes de-projected in the ^-direction to show how they would appear on the sky (i.e. Fig. 1131 The axes are plotted in 
co-moving distance (Mpc) so the boxes have sides of RT00 Mpc (except the central figures which have a larger extent in the z direction). 



4-7.1 Press- Schechter approach 



If the bias factor were ~1, the super-structures could be 
treated as non-linear collapsing objects and we could then 
use a Press-Schechter approach. Sec. 14.61 shows the radio 
galaxies in both super-structures trace an approximately 
spherical distribution which is an assumption in this ap- 
proach. To apply the Press-Schechter appr oximation, we fol- 
lowed the method of ISteidel et al.l il998T) . Assuming some 
value for the bias factor b, we calculated the mass of each 
super-structure given by M =< p > V(l + <5 m ) where < p > 
is the mean density of the Universe, V is the co-moving 
volume and <5 m is the mass overdensity. We find a mass of 
M ~ 7 x 1O 17 M0 for each super-structure. We then calcu- 
lated the probability of observing a peak this high in the sur- 
vey volume in a ACD M Universe using the Press-Schechter 
approximation given in lSteidel et alJ {1998). For a bias fac- 
tor, 6=1, we obtained a probability of < 10 -11 . We conclude 
that, like other radio galaxies, the TONS08 radio galaxies 
are not unbiased tracers of mass. 



4-7.2 Quasi-linear structure formation theory 

As we increase the bias factor, the radio galaxy overden- 
sity corresponds to a smaller mass overdensity. Structures 
go non-linear at overdensities of approximately 1.7. There- 
fore, as the bias increases, the super-structures are still in 
the quasi-linear regime (albeit only just) and the Press- 
Schechter method is not applicable. As discussed above, 
we have good reason to expect the bias factor to be 
greater than 1 as locally , it has been measured to be ~1.65 
^Peacock fc Doddslll994) . We therefore adopt the following 
method. 

We model the height of the cosmic power spectrum 
at the relevant scale. We consider the dimensionless power 
spectrum, A 2 (fc) as this gives us the dimensionless spectrum 
of density fluctuations per In fc. As we don't know the full 
extent of the structures, we calculate this for the scale 50 
Mpc (radius) (fc=0.126). The normalisation of the power 
spectrum is fairly well known down to wave numbers of ap- 
proximately 0.028 dPeacock fc Dodddll994h . 

We calculate the power spectrum at redshift z and for 
a wavenumber fc using 
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(11) 



where a = 1/(1 + z), G z is the growth factor at redshift z 
feisenstein fc Hulll998f) . is the transfer function as given 
in iBond fc Efstathioul (|l984), n is the index of the power 
spectrum and q = k/T. The effective wavenumber k B g is 
defined as the wavenumber at which erg measures power, 



k eS /h = 0.172 + 0.011 



In 



(0.34) 



(12) 



JPeacockll200ll) . This then defines q$ = k e g /T and a transfer 
function t$ calculated at fc e ff. We then estimate the first 
three moments of the variance 



a 2 m = / A 2 (k)\ f k iV^dk, 



(13) 



where | /fc | is the Fourier transform of a 3-D top-hat fil- 
ter function in which we set the radius to 50 Mpc and m is 
the index of the moment. For the first moment, we obtain 
a value of 00 = 0.243. This is the rms density variation in 
50 Mpc spheres. We can check this value is correct by the 
following rough calculation. Adopting the canonical form of 
the correlation function for galaxies, £ = (r/ro) -1 ' 8 , and as- 
suming <7q w £, we can calculate a^ = o% h -i x (50fe/8) -1 ' 8 , 
where cr gh -i = 0.94. We obtain the value 050 = 0.249 which 
agrees well with the above ca lculation. 

The spectral parameters feardeen et aljll986l) are then 
calculated 



v= — ,7 = — — ,7?* = V3— , (14) 

(To CT QO~2 CT2 

5m is estimated as before using 8 m —5 sa ,\/b. The rarer the den- 
sity peak, the more it will have significant departures from 
a simple Gaussian random field due to non-linear growth 
by self-gravitation. W e use a log-normal (LN) random field 
<Coles fc Jones! Il99ll) to make some attempt at correcting 
the mass overdensity, 8 m to a non-linear corrected mass over- 
density Snl 



<5 n i = log(l + <5 m ) + ■ 
where 

"nl = log (O"o + 1). 



(15) 



(16) 



The probability of intercepting a structure of height 
v or greater was calculated by multiplying the co-moving 
volume of TONS08 (1.5 x 10 7 Mpc 3 ) by the cumulative 
number density of peaks higher tha n v. We used the fit- 
ting formulae of lBardeen et all il98rt) (equations 4.3,4.4,4.5 
and 4.11a). The results are shown in Fig. 1151 As the 
results are very dependent on the normalisation of the 
power spectrum, we performed the analys is for g s> ,-i= 0.75 
from combined 2dF/CMB measurements JSeliak 2002f) and 
aft,, -i=0.94 from weak gravitatio nal lensing measurements 
jRefreeier. Rhodes fc Grothl2 0021. The results are shown in 
Fig. [Inland Table. H 

Without taking into account non-linear effects, if we as- 
sume a relatively high value for cr sh -i (0.94) and that the 
radio galaxy bias parameter is the same as for smaller scales 
in the local Universe (1.65), we have found an extremely 
rare (8cr) fluctuation (the 2=0.27 super-structure). If we in- 
clude a non-linear correction, the effective mass overdensity 



is reduced and corresponds to a much less extreme rare fluc- 
tuation of 3.9cr. After non-linear corrections, we calculate 
that we would only expect 0.05 such super-structures in our 
sam ple volume. These co rrections are therefore very impor- 
tant. ISteidel et aD <199S^ use a different method for making 
non-linear corrections jBernardeaul|l994h which give similar 
results. 



4-7.3 Redshift space distortions 

In calculating the radio galaxy overdensity, we may need 
to account for redshift space distortions. As 8p/p is close to 
the critical value, we have until now assumed that the super- 
structure is just breaking away from the Hubble expansion. 
However, if it is more evolved and is collapsing, then the 
super-structure will appear more dense than it actually is. 

If we assume a maximum peculiar velocity v pcc towards 
the centre of the superstructure, we can work out the actual 
redshift of each radio galaxy and hence how many radio 
galaxies are actually outside the 50 Mpc radius sphere and 
calculate the new radio galaxy overdensity. 

The maximum radial infall velocity i> pcc around a spher- 
ical perturbation of radius R and interior average mass over- 
density (5 m is given by 



;H Rf5 n 



(17) 



where fw SI™ l|Lahav et al.lll99ll ). 

The problem in calculating this maximum infall velocity 
is that as we increase the value of the infall velocity, fewer 
radio galaxies fall within the 50 Mpc radius sphere. This in 
turn decreases the mass overdensity which will decrease the 
peculiar velocity (due to Equation 17) . We therefore choose 
to find a compromise value for the redshift-space distortion 
corrected radio galaxy overdensity 5™f by calculating the 
radio galaxy overdensity given different v pcc for both the 
data and theory (from equation 17 and assuming a radio 
galaxy bias factor of 6=1.65). We find a "corrected" radio 
galaxy overdensity of (5g°[=1.35 (10 radio galaxies) and 1.84 
(8 radio galaxies) for the 2=0.27 and 0.35 super-structures 
respectively (see Fig. 11611 . Of course, real collapsing systems 
will be more complex and hence have more complex velocity 
fields, so this is only a rough calculation. 

Using these corrected radio galaxy overdensities (and 
assuming a radio galaxy bias factor 6=1.65), the spectral 
parameter u—2.84 and 3.58, the non-linear corrected mass 
overdensity 5 n i=0.72 and 0.87, and the expected number of 
peaks of height v or greater in the survey volume =0.87 and 
0.12 for the 2=0.27 and 0.35 super-structures respectively. 
Therefore, if these super-structures are collapsing, this has 
a huge effect on the results and it is not that surprising that 
we detect them in the TONS08 survey volume. 



5 DISCUSSION 

5.1 The 2=0.27 super-structure 

5.1.1 The geometry of the z=0. 21 super- structure 

The spatial extent of the 0=0.27 super-structure in TONS08 
is w 80 x 100 x 100 Mpc 3 (see Sec. EHJ. A turn over in 
the matter power spectrum at these scales is still disputed 
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No. density No. density 



Figure 15. The probability of intercepting a super-structure of radius 50 Mpc and height associated with the galaxy overdensity of 
the ,2=0.27 super-structure (left) and z=0.35 super-structure (right) against the radio galaxy bias factor. Og^-i is assumed to be 0.94. 
Overplotted are the 68 per cent confidence limits of the galaxy overdensities which gives us an indication of the effects of Poisson sampling 
fluctuations. A local value for the radio galaxy bias parameter b=1.65±0.15 is plotted by horizontal lines with the regions within the 
error shaded. 
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Table 4. Table showing the mass overdensity 5 m , the non-linear corrected mass overdensity <5 n i, the mass associated with this overdensity, 
the spectral parameter v (the height of the field in units of the rms) and the expected number of peaks of height v or greater in the 
survey volume for different values of radio galaxy bias and a sh -i for the two super-structures. 



{Miller fc Batuskill200lh . mostly due to problems with win- 
dow functions and sample incompletenesses. Recent results 
from the full 2dfQSO survey suggests no t urnover in the 
powe r spectrum below scales of 210 Mpc ijOutram et alJ 
2003). However, it is clear that our radio-selected super- 
structures are comp arable to the largest structures found 
to date. For example iBatuski et alJ jlQQOT) use optical clus- 
ters to find four super-clusters with a maximum extent of 
~80 Mpc. These structures, as well as most structures found 
to date, do not extend as far in all directions as TONS08 



(although of course radio galaxies trace the TONS08 super- 
structures so sparsely, it is difficult for us to make firm state- 
ments) . 

Sec. 14.61 makes an attempt to quantify the 3-D distri- 
bution of radio galaxies within the super-structures. The 
z=0.27 super-structure appears more elongated roughly in 
the z-direction. If we assume that the super-structure is in- 
falling, then redshift space distortions (se e Sec. 14. 7. 311 will 
make this an even more pronounced effect. IZeldovichl jl970tl 
argued that pancakes are the first structures formed by grav- 
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Figure 16. The redshift-space distortion corrected radio galaxy ovcrdcnsity for different values of maximum peculiar velocity for both 
the 2=0.27 (left) and 2=0.35 super-structures. The results are shown for the data (solid line) and theory (dotted line). We chose a 
compromise radio galaxy overdensity by determining where these curves cross. 



itational clustering and are the dominant features arising 
from the first stages of non-linear gravitational clustering. 
They flatten out into sheets and then break up into filaments 
or form directly into sheets. The 2=0.27 super-structure 
could therefore plausibly be elongated in a particular di- 
rection because of the way it is collapsing. Alternatively, it 
is possible that the structure is bigger in all three dimen- 
sions, but the angular extent of the TONS08 survey means 
that we don't see all of it. If true, then this structure would 
become truly remarkable as it would correspond to a very 
high-sigma fluctuation 2 . It is also plausible that what we are 
seeing is actually two super-structures which are in-falling 
into each other along our line of sight. Such double super- 
structures may not be that unusual for the same reason, 
nam ely bias, that mergi ng clusters are relatively common 
(e.g. iBardelli et alll2002l) . 

We will now discuss ways of reconciling the low prob- 
ability of TONS08 intercepting the 2=0.27 super-structure 
with the fact that we have detected it. It is possible, ar- 
guably probable, that two or more of these explanations are 
important. 

5.1.2 The value of a sfl -i 

These results are very dependent on the assumed normali- 
sation of the power spectrum, a 8h -i. If we assume a lower 
value for cr gh -i as preferred by combined 2dF/CMB mea- 
surements, then the expected number of super-structures in 
TONS08 decreases to 9xl0~ 4 . There is still much contro- 
versy over the value of a sh -i. Recent work suggests that 
the higher values of cr 8h -i obtained from weak lensing tech- 
niques may be due to systematic errors and that the values 

2 We have therefore conducted a new, wider angle, higher radio 
flux density survey (TOOT08w: Table. HI to further constrain the 
dimensions of the super-structure (Brand et al. in prep.). 



are converging to lower values in line with studies of X-ray 
clusters and combined 2dF/CMB estimates. However, these 
results all assume fi m = 0.3. Assuming a slightly lower value 
for f2 m would increase o~ 8h -i to nearly 1. In all of this anal- 
ysis, our results actually depend on o~ sh -i x b. We there- 
fore favour (j gh - 1=0.94 in this paper because it is closest to 
o~ sh -i x b for optical galaxies and is a much better deter- 
mined quantity than that of the dark matter. 

5.1.3 A different radio galaxy bias factor? 

Another explanation of our results is that radio galaxies 
are extremely biased tracers of mass in the 2=0.27 super- 
structure (i.e. more biased than local radio galaxies). This 
could be due to redshift evolution of the bias factor or to 
different mechanisms being at work on larger scales, or a 
combination of the two. 

An evolution of bias with redshift could be linked to the 
well known evolut ion in the space density of radio galaxies 
with redshift (e.g. IWillott et alJl20Q]]) . ICroom et all fcOOll) 
look at the evolution of QSO clustering as a function of 
redshift in the 2dF QSO redshift survey. They find that 
QSOs have similar clustering properties to local galaxies 
at all redshifts in their survey (z^3.0). Because the clus- 
tering of dark matter is thought to decrease at higher red- 
shifts (because fluctuations grow under gravity), this implies 
an increase in the bias factor of optically-selected QSOs at 
higher redshifts. The si milar clustering properties o f radio 
galaxi es at low r edshifts aPeacock fc Nicholsonlll99ll) and at 
2 ~1 (Blake & Wall 2002) suggests similar results for radio 
galaxies. If the correlation length of radio galaxies is simi- 
lar at high redshifts to that found locally, there must be an 
evolution in bias. Assuming a linear evolution of bias with 
redshift of the form b(z) oc (1 + z) n , a value of n=0.4 is re- 
quired for stable clustering in proper coordinates. Taking the 
local radio galaxies bias factor to be 6(0) = 1.65, this implies 
a radio galaxy bias factor 6(0.27) = 1.8 at 2=0.27. This in- 
crease in the bias factor with redshift is fairly small at these 
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relatively low redshifts and is also very model dependent. 
However, because the probability of TONS08 intercepting 
a super-structure is such a steep function of bias, only a 
small increase in bias may be enough to mean that the ra- 
dio galaxies are tracing a less rare fluctuation in the evolved 
power spectrum that we would be more likely to intercept. 
We note here however that these calculations are based on 
data within a much larger survey area. A larger redshift evo- 
lution of bias may be possible if the triggering mechanism 
that causes a higher bias is peculiar to super-structures. We 
now discuss the possibility of a different radio galaxy bias 
for large scales. First, the 2=0.27 super-structure looks to 
be an unusually dense region of sky. It is likely to be rich in 
collapsed substr ucture, e.g. rich clusters of galaxies (see also 
| Hill_el_aD |20Q2|) . The bias factor for rich cl usters is ~3.1 
JPeacock fc Dodddll99l IPlionis et all |2000T) . This is pre- 
cisely the value required to make the probability of seeing 
a super-structure in TONS08 ~ 1 (Fig. I15fl which would be 
the case if the radio galaxies are associated with rich clus- 
ters. Enhanced radio galaxy triggering due to group-group 
mergers (e.g. cluster forma tion) should also be considered 
( Sim pson fc Raw lings 2002). Large surveys of radio galax- 
ies would include few of these rare systems (i.e. clusters or 
forming clusters) and on average give lower valu es of bias 
JPeacock fc Nicholsonl Il99ll iBlake fc Wallll2002l) i.e. more 
appropriate to typical radio galaxy environments. 

The fact that the system may be collapsing may 
also be relevant. At recent epochs, structures the size 
of rich clusters are going non-linear. At least the cen- 
tral parts of s uper-structures have reached t heir turn- 
around radius jGramann fc Suhhonenko| 12002ft so per- 
haps this effect is peculiar to collapsing systems. Re- 
cent observations of large-scale structure in the 2dFGRS 
iNorber ge^lJ|2002j), N-body simulations and cluster cat- 
alogu es feolokoko. Basilakos fc Plionisll2002l . lEinasto et alJ 
1997) reveal that the matter in the Universe follows a fila- 
mentary structure. The largest structures lie at the nodes 
which join filaments and material falls into these along 
the filaments. Matching radio galaxies with 2dF galaxies 
reveals that radio galaxies are often on the edge of high 
density regi ons (Brand et al. in prep.) This is a lso seen 
for quasars (ISochting. Clowes fc Campusanoll2002l) . In the 
TONS08 z=0.27 super-structure, only 2 of the 4 7CRS (the 
radio-brightest objects) are in the central region defined by a 
50 Mpc radius sphere. Perhaps radio activity is triggered by 
the high velocities as the galaxies fall into the centre of the 
super-structure. Indeed, there are suggestions that nearby 
radio galaxies such as those in the Pers eus cluster may hav e 
been triggered by their high velocities JPedlar et al.l ll990). 
Note, also, that some degree of redshift space distortion is 
inevitable in a collapsing system which will inevitably and 
spuriously enhance the apparent overdensity of the super- 
structure (see Sec. 14. 7.31 . 

5.1.4 Poisson sampling fluctuations 

Fig. 1151 also shows the 68 per cent confidence limits of 
the galaxy overdensities on the probability of TONS08 in- 
tercepting a super-structure as a function of bias. Poisson 
sampling fluctuations could mean that there happen to be 
more or fewer radio galaxies than we would expect on av- 
erage for a fluctuation of a given size. This could modify 



our results slightly in either direction but not enough to 
explain our results. Fig. 1101 shows that this is a more se- 
vere problem for sparser samplers of the mass overdensity, 
such a s radio gala xy surveys with higher radio flux density 
limits llLacvll2000l) or radio galaxy surveys at low redshifts 
at which the space density of radio galaxies is much lower 
JPeacock fc Nicholsonlll99lh . 

5.1.5 A representative volume? 

In the previous analysis, we have assumed that the TONS08 
survey area is completely random. Although the original 
three 7CRS fields were essentially random in their position 
on the sky, our reason for choosing the TONS08 survey area 
was to follow up an overdensity already noticed in the 7C-II 
region. This obviously makes it more likely that we found 
an overdense region but by a factor that is hard to quantify 
until we have complete large-scale structure measurements 
in other TONS regions. 

5.2 The z=0.27 and z=0.35 super-structures 

If we again assume a ah -i=0.94 and that the radio galaxy 
bias parameter 6=1.65, Fig. ll5l forces us to conclude that the 
0=0.35 super-structure is an even rarer (5.1a) fluctuation 
of which we would only expect ~3.3E-4 in our sample vol- 
ume. Assuming that it is real, the same arguments as made 
in Sec. 15. II apply to the 2=0.35 super-structure. Because it 
corresponds to a rarer fluctuation, and because there are 
two super-structures in the volume, many of the arguments 
are not strong enough. For example, a simple evolution of 
bias with redshift would have to be fairly extreme and would 
contradict other studies (e.g. lBlake fc Wallll2002T) . 

An analysis of angular clustering in TONS08 (Sec. 14.51 
shows a possible edge in the distribution of radio galaxies 
within the 2=0.35 super-structure. Its significance (Sec. 14.21 
is determined for the whole angular extent of TONS08. This 
will not take into account the fact the the 2=0.35 super- 
structure is probably smaller than this and that number 
of radio galaxies expected in the same volume will also be 
smaller. Fig. llOl shows that it is fairly unlikely that the over- 
density is small enough to make the collection of radio galax- 
ies not significant. For the following discussion, we therefore 
assume the super-structure to be real. 

Perhaps, the two fluctuations are linked. They are sep- 
arated by about 310 Mpc and the power spectrum may sti ll 
have power on these large scales (IMiller fc Batuskil l2001). 
The rarer the fluctuation, the more highly clustered it will 
be. Perhaps we observe these two rare fluctuations because 
they are sitting on top of an even rarer larger-scale fluc- 
tuation. The two-point correlation function is defined as 
the excess over Poisson of the joint probability of finding 
objects in two volume elements which are sepa rated by a 
distance r. Using equation 3 from Kaiser (1984), we deter- 
mine that the correlation function for the super-structures 
is boosted by a factor of «300 over the typical value at this 
separation. Although boosted by a large amount, the cor- 
relation function of mass at these separations is very small 
(~0.001). The correlation function will only be boosted to 
a modest 0.5. Therefore, the probability of finding a second 
super-structure is 1.5 times higher than it would otherwise 
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be (half of all super-structures separated by 300 Mpc are 
genuinely clustered and half are there by chance). This is a 
fairly negligible effect compared to the other uncertainties 
so we can effectively just multiply together the probabili- 
ties of finding each of the super-structures to obtain a joint 
probability of ~ 10 _J . As we discussed in Sec. 15.11 we may 
have targeted an especially dense region of the Universe be- 
cause we h ave followed up a stru cture already discovered in 
the 7CRS jRawlines et al]l2003ft . If this is the case, we may 
only be able to include the newly discovered 2=0.35 super- 
structure in our analysis, bringing the joint probability down 
to ~ 10" 4 . 

We conclude that an increase in bias on the large scales 
of the super-structures and/or redshift space distortions are 
probably needed to explain finding two super-structures in 
the TONS08 region. 

5.3 Finding similar super-structures in other 
surveys 

iGraham. Clowes fc Campusanol 1^995) combine the results 
of six quasar surveys to look for super-structures and find 
evidence for three structures. These are sm aller and more 
elonga ted than the structures we have found. IWillieer et alJ 
(2002) find evidence for a huge structure ~ 70 x 140 x 140 
Mpc 3 as traced by 18 quasars at redshift 2=1.3. This is 
the largest st ructure at high red shifts yet found. Both 
the 2dFGRS JColless et alJ boOlft and 2dFQSO surveys 
JCroom et al]l200ll) should be able to find super-structures 
especially, if radio emission is a particularly biased tracer 
of mass, when combined with radio surveys such as NVSS 
or FIRST. The 2dFQSO survey has some evidence for clus- 
tering on large scales. However, they have not yet reported 
evidence for any super-structures. Radio loud quasars have 
been studied but only approx 2 per cent of optical and ra- 
dio IDs match up meaning that the space density of radio 
loud quasars may not be sufficient at z ~0.5 to trace super- 
structures. Similarly, the 2dFGRS team have not yet re- 
ported any newly discovered super-structures. Possible rea- 
sons include the fact that matching algorithms may cause 
incompletenesses, sky incompletenesses may prevent stud- 
ies of the largest scales or it may be simply that noone has 
looked properly yet. 

Sec. 14. 3. T1 discusses the discovery of 2 redshift peaks in 
the lLacvl i2000t) survey. The corresponding mass overden- 
sities (assuming 6=1.65) are so large that even though the 
survey volume is w 7-times larger (i.e. 7-times more inde- 
pendent 50 Mpc radius spheres), the probability of there 
being two such super-structures in the survey are tiny (~ 
1E-9 for the 2=0.28 redshift peak). We should be cautious 
of this result for the following reasons. Because of the small 
number statistics, the errors are larger than for TONS08 and 
there is a good chance that the real overdensities are smaller 
(Fig. I1UI and hence correspond to less rare primordial fluc- 
tuations. The log-nor mal approximation fo r the non-linear 
corrections we apply JColes fc Joneslll99ll) are also known 
to break down for the highest peaks so the non-linear cor- 
rected overdensities may be inaccurate. Even if these super- 
structures are only as overdense as the TONS08 superstruc- 
tures, the probability of intercepting two such structures is 
still fairly low. lBennl <ll997f) finds a group of four radio galax- 
ies at 2 «0.31 spanning ~20 Mpc on the sky in a similar 



radio galaxy redshift survey. It seems that these huge struc- 
tures in other surveys may not have looked very significant 
because the survey area is less well matched than TONS08 to 
super-structure scales and/or because of the effects of Pois- 
son sampling fluctuations when the average space density of 
the radio galaxies are low. 

If an increase in bias is not due to a redshift evolution 
but instead due to some specia l triggering mechanism in col- 
lapsing super-structures then iPeacock fc Nicholson! Jl99lfl 
should have detected similar super-structures. After all, even 
in a ACDM Universe, there are expected to be more super- 
structures of a given mass collapsing per unit volume at 
lower redshift. With a redshift limit of z ~0.1 and sky area 
of 9.3 sr, their survey volume is naively of order 15-times 
larger than the volume probed by TONS08. They find 9 
structures with 4 or more radio galaxi es within a 57 Mpc ra- 
dius s phere. The flux density limit of IPeacock fc Nicholson! 
Jl99lf) is Si. 4 -500 mjy with a redshift limit of z wO.l. If 
we compare this to TONS08, we find that the radio lumi- 
nosity of the faintest radio galaxy detectable at the redshift 
limit at z « 0.1 is w 5x that of TONS08. By calculating 
the model redshift di s tribut ion for the survey limits of the 
IPeacock fc Nicholson! Jl99lJ) sample, we predict that there 
should be a total of 312 radio galaxies at z <0.1 (this com- 
pares well with the total number of 310 in the sample). With 
a co-moving volume of 2.25 xlO 8 Mpc 3 , we expect 290 in- 
dependent spheres of radius 57 Mpc. We would therefore 
expect a mea n number of rad i o gal axy per 57 Mpc 

radius sphere. IPeacock fc Nicholson! Jl99lT) only look for a 
minimum of 4 objects in each 50 Mpc radius sphere. As 
only 3.3 radio galaxies are needed to produce an overden- 
sity of 2, the small number statistics involved at low red- 
shifts ensure that many (m 60 per cent of) overdensities of 
this size will not be detected. Using the method in Sec. 14.71 
and assuming the canonical bias of 1.65, we would expect 
IPeacock fc Nicholson! dl99lF) to have found ~ 0.8 of these 
in their survey (assuming a redshift o f ~ 2=0.05). Because 
they l ook for a minimum of 4 objects, IPeacock fc Nicholson! 
( 1991) should only have been able to detect overdensities 
of 2.6 or more (neglecting Poisson uncertainties). Again as- 
suming the canonical bias of 1.65, we would expect them to 
have found only 1.7E-2 of these in their survey volume. If 
we increase the bias to &=4, our model predicts that they 
should find 11 which is in good agreement with what they 
find. 

The fact that similar structures appear in other surveys 
at both similar and lower redshift, reinforces the suggestion 
that we are seeing some mechanism at work which boosts 
the bias of the radio galaxies within super-structures. 

5.4 Alternative structure formation theories 

The currently accepted theory of structure formation is that 
of the freezing in of quantum fluctuations of a scalar field 
during an inflationary period. An alternative theory that 
predicts larger numbers of rare fluctuations than the ac- 
cepted inflationary theory is that large-scale structure is 
seeded by topological defects which are formed naturally 
during a symmetry br eaking phase transition in the early 
Universe iTuruklfl996h . Although recent measurements of 
the po sition and amplitu de of the acoustic peaks in the 
CMB dBenoit et alJ 2003) rule out cosmic defects as the 
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dominant structure formation mechanism, it is possible that 
they provide a non -negligible contribution at large scales 
l)Sakellariadod r2002V This would boost the number density 
of super-structures. 



6 CONCLUDING REMARKS 

In this paper we have presented evidence of two significant 
super-structures in the TONS08 radio galaxy redshift survey 
as traced by radio galaxies in 50 Mpc radius spheres. From 
the implied overdensities of m 3 on these large scales, we 
calculate the mass of these super-structures to be « 5 x 
1O 17 M0. The super-structure at 2=0.27 (traced by at least 
13 radio galaxies) w as tentatively identified in the 7CRS 
jRawlines et all 2003 ) , traced by the brighter radio galaxies. 
The super-structure at z=0.35 (traced by at least 12 radio 
galaxies) is newly discovered. 

Assuming <r gh -i=0.94 in a ACDM Universe and a radio 
galaxy bias factor of 1.65 these super-structures correspond 
to the evolved counterparts of ~ 4 and 5cr fluctuations in 
the primordial density field. The expected number of such 
fluctuations in a volume equivalent to that of the TONS08 
survey is ~ 0.05 and ~ 3 x 10~ 4 . 

The probability of seeing two such fluctuations in 
TONS08 is strongly dependent on how representative the 
survey region is of the Universe. If, by targeting a structure 
discovered in the 7CRS, we have focused on a particularly 
unusual part of the Universe, then the probability of seeing 
two super-structures is roughly equal to the probability of 
finding only the z=0.35 super-structure, namely ~ 10 -4 . If, 
however, the TONS08 target region turns out to be a rea- 
sonably representative one, then given our assumptions, the 
probability of seeing two super-structures is ~ 10 . The 
presence o f similar radio galaxy overdensities llLacvll2000l 
lBenrJl997ft at similar redshifts suggests that these overden- 
sities are actually fairly common. 

From Fig. 1151 it is clear that the probability of detect- 
ing a radio-selected super-structure is a very strong function 
of the assumed radio galaxy bias. If we assume evolution of 
bias with redshift of the form b(z) oc (1 + z)", we would need 
n ~ 2 to increase the probability of intercepting a super- 
structure such as the z=0.35 super-structure to a reason- 
able value. Previou s studies suggest that this is an unfeas i- 
bly large evolution JCroom et all200ltlBlake fc Wallll2002>) 
Simila r structures are also found by iPeacock fc Nicholson! 
Jl99ll) at low redshifts, again suggesting that this cannot be 
a purely evolutionary effect. 

Alternatively, an increase in bias could be due to some 
special triggering mechanism in collapsing super-structures. 
This seems the most promising way of reconciling the low 
probability of TONS08 intercepting both super-structures 
with the fact that we have detected them. There are sev- 
eral plausible explanations for this. Perhaps collapsed sub- 
structure within the super-structure, i.e. rich clusters, hosts 
the radio galaxies and the high bias simply reflects that of 
the rare rich clusters. However, sim ilar super-structures have 
been found with radio quiet QSOs (Wil lieer et alJl2002fl and 
these tend to be found in less rich environments. Perhaps 
there is an enhanced rate of group-group mergers within 
the super-structure which is enhancing radio galaxy trigger- 
ing. Perhaps as the super-structure goes non-linear, redshift 



space distortions and/or enhanced triggering due to high 
velocities boost the bias. 

We have just completed radio galaxy redshift surveys 
in other reg ions. These, along with radio galaxies selected 
from SDSS JStoughton et al.ll2002l) . should determine how 
rare the TONS08 super-structures really are and, assum- 
ing structure formation models are correct, accurately deter- 
mine the radio galaxy bias within super-structure regions. 
To determine the origin of this enhanced bias, we have begun 
a K-band (UKIRT) imaging survey of radio galaxies within 
the super-structures. By comparing structural parameters, 
merging rates and clustering environments with radio galax- 
ies of similar optical magnitudes and radio flux densities out 
of super-structures, it should be possible to discriminate be- 
tween the above possibilities. 

We conclude that despite the seemingly remarkable dis- 
covery of two rare, huge and massive super-structures in 
the fairly small volume of the TONS08 survey, this proba- 
bly does not present any serious challenge to the standard 
ACDM model and the inflationary (i.e. Gaussian fluctua- 
tion) paradigm for structure formation. 
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08 
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24 
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33 
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19 
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TONS08.880 


NOT 


2001 
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08 


24 


39, 
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26 
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2.7 
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08 


25 
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,13 


27 
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26 


,84 
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08 


25 


54. 


,94 
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04 
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.32 
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2.7 
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Mar 14 


08 


26 


37. 


.06 


26 


22 


58 


.86 


TONS08.1020 


WHT/2.7 
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Apr 17/Nov 7 


08 


27 


07. 
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28 


10 


22 


.06 


TONS08_1031 








08 
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.40 


29 
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.64 
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NOT 
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08 
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.84 
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NOT 
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08 
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WHT 
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08 
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28 


45 


18 
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08 


28 
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,54 


27 


47 


02 


.14 


TONS08.1099 


2.7 
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08 


28 
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,33 


24 


20 
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NOT 
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Oct 31 


08 


28 


42. 


,17 


24 


30 
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.74 
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NOT 
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08 
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29 


17 


21 
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NOT 
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Oct 31 


08 


28 
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,80 


24 


42 


39 
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TONS08.1124 








08 


28 
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,07 


24 


33 


25 


,34 


TONS08.1140 
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Oct 31 


08 


29 
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26 


00 


01 
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Z ' 
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7.8 


7.82 


18 


.45 


21 


.78 


0.366 





.3654 





.3662 




6.7 


5.94 


18 


,51 


20 


.69 


0.260 





,1985 





.1994 


[OII]/[OIII] lines 


25.7 


23.74 


19 


,19 


20 


.26 


0.453 





,4522 





.4536 




14.4 


5.34 


18 


.37 


21 


.66 


0.337 





,3374 





.3364 




9.7 


10.74 


19 


.21 


21 


.48 


0.393 





,3927 





.3917 




13.0 


10.15 


17 


.51 


20 


.44 


0.280 





,2798 





.2800 




4.5 


3.52 


18 


.99 


21 


.78 


0.347 





,3458 





.3451 




5.3 


2.27 


18 


,70 


21 


.38 


0.357 





,3579 





.3578 




7.5 


5.91 


19 


.68 


21 


.78 


0.442 





,4405 





.4418 




4.8 


5.51 


16 


.80 


19 


.51 


0.315 





,3150 





.3188 




3.4 


2.92 


17 


.83 


21 


.48 


0.260 





,1580 





.1650 


[OII]/[OIII] lines 


8.3 


8.06 


19 


.38 


21 


.78 


0.363 





,3618 





.3619 




144.7 


133.57 


15 


.42 


17 


.14 


0.268 
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5.6 


6.42 


14 


.28 


18 


.34 


0.174 





.1740 





.1750 




22.7 


6.72 


18 


,28 


21 


.78 


0.281 





,2810 





,2815 




21.3 


9.14 


19 


,63 


21 


.48 


0.390 





.3910 





.3907 




6.6 


6.06 


19 


,44 


21 


.48 


0.448 





,4480 





.4478 




70.3 


3.41 


19 


,06 


70 


.34 


0.347 
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20.1 


15.47 


19 


,05 


21 


.78 


0.345 





,3439 





.3432 




22.5 


10.26 


16 


,52 


19 


.29 


0.136 





.1371 





.1365 




3.8 


3.67 


18 


.68 


21 


.48 


0.397 





.3970 





.3974 




3.7 


3.38 


15 


.05 


17 


.32 


0.078 





.0816 





.0807 


Ha line 


8.4 


7.35 


18 


.17 


20 


.44 


0.251 





.2502 





.2496 




15.9 


11.35 


19 


.60 


21 


.78 


0.401 





.3997 





.3986 




6.5 


5.18 


19 


.54 


21 


.48 


0.357 





.3564 





.3577 




5.7 


8.10 


16 


,07 


18 


.81 


0.103 





.1031 





.1031 




5.7 


2.28 


19 


,35 


21 


.78 


0.370 


0, 
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17 
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.83 
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2043.1 
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19 


,15 


21 


.78 
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10.4 


6.83 


18 


.67 


21 


.78 
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6.43 


18 


.15 


21 


.48 


0.381 





,3815 
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17 


,70 


21 
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15.71 
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,31 


21 


.78 
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,3688 





.3692 
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1.23 


15 


,80 


18 


.19 


0.136 





.1357 
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10.7 


3.60 


18 


,49 


21 


,48 


0.361 





.3611 





.3599 




9.0 


8.80 


19 


,83 


21 


.78 


0.239 





.2384 





.2381 




5.1 


2.12 


18 


,33 


20 


.50 


0.268 





.2679 





.2680 




123.4 


8.82 


17 


,76 


37 


.35 


0.243 
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9.07 
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